Critical exponents at the ferromagnetic transition in 
tetrakis(diethylamino)ethylene-C6o (TDAE-Cgo) 
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Critical exponents at the ferromagnetic transition were measured for the first time in an or- 
ganic ferromagnetic material tetrakis(dimethylamino)ethylene fullerene[60] (TDAE-Ceo)- From a 
complete magnetization-temperature-field data set near Tc = 16.1 ± 0.05, we determine the suscep- 
tibility and magnetization critical exponents 7 = 1.22 ± 0.02 and {3 = 0.75 ± 0.03 respectively, and 
the field vs. magnetization exponent at Tc oi 5 = 2.28 ± 0.14. Hyperscaling is found to be violated 
by f2 = d' — d ~ — 1/4, suggesting that the onset of ferromagnetism can be related to percolation of 
a particular contact configuration of Ceo molecular orientations. 



Molecular ferromagnetism, particularly when only 
electrons in p orbitals are involved in the magnetic in- 
teractions, is a relatively newly discovered phenomenon. 
Of the few compounds discovered so far which display 
signatures of proper ferromagnetic (FM) behavior , 
the most studied has been TDAE-Cgo, which - by virtue 
of its relatively simple synthesis (at least in powder form) 
and high Curie temperature of Tc =16 K - has been in- 
vestigated by many groups. However, the magnetic prop- 
erties of this material are not straightforward, and mea- 
surements on powder samples have lead to apparently 
conflicting proposals regarding its low-temperature state 
ranging from superparamagnetism |^ to spin-glass ^] as 
well as ferromagnetism. More recently, low-field electron 
paramagnetic resonance 1^ and susceptibility measure- 
ments on high-quality single crystals showed more con- 
clusively that the material - if properly chemically and 
thermally prepared |^ - display clear signs of a transition 
to a ferromagnetic state at around 16 K, in agreement 
with the original suggestion of its discoverers |0] . 

The origin of the ferromagnetic exchange interaction 
between Ceo molecules in TDAE-Cgo has been studied 
theoretically by a number of groups Recently it 

was discovered Q that in the FM a-phase of TDAE- 
Ceo two different orientations of Ceo-molecules may oc- 
cur at low temperatures (labeled / and //). These ori- 
entations can lead to different contact Ceo configurations 
along the direction of closest contact (c-axis) ||], which 
profoundly affect the exchange interaction along c-axis 
0. In the non-magnetic a' phase, the 6-6 double bond 
(nearly) faces the center of the hexagon on the neighbor- 
ing molecule, whereas in the FM phase, a number of ad- 
ditional different mutual orientations are possible. But, 
among these, the alternating /-// contact configuration 
- in which the double bond on one molecule approxi- 
mately faces the center of the pentagon of its neighbor - 
was shown to be dominant in the ferromagnetic state. 

To confirm the observation of a proper FM state and 



investigate the associated critical behavior, we report 
here the first measurements of the critical exponents as- 
sociated with the ferromagnetic transition of TDAE-Ceo- 
The measurements, which are also the first for any or- 
ganic system, are found to give remarkably self-consistent 
values of the critical exponents in agreement with behav- 
ior expected for a ferromagnetic transition in a system 
with a certain degree of disorder. They confirm the es- 
sentially FM behavior, and also give important insight 
into the interactions responsible for the ferromagnetism 
in these unusual materials. 

Single crystals of TDAE-Cgo where grown by a diffu- 
sion method as described in ||^. For magnetic measure- 
ments, a number of crystals from different growth batches 
were sealed into quartz tubes under helium. Since the 
"as grown" crystals of TDAE-Ceo are in their a' mod- 
ification which shows no low-temperature ferromagnetic 
transition |p^ , they were annealed in order to transform 
them into the ferromagnetic a modification. The an- 
nealing was carefully done through several intervals of 
1 hour at 70° C, each of them followed by a measure- 
ment of the low-temperature magnetic properties. The 
annealing procedure was stopped at the point were the 
low-temperature saturation magnetization of the crystals 
reached it's maximum, i.e., when the whole sample had 
transformed into the a modification. Magnetic measure- 
ments were performed with a Quantum Design MPMS 
SQUID magnetometer which enables a temperature sta- 
bility better than 10 mK and a measurement of magneti- 
zation with a relative error of less than 0.1% and since the 
magnetic response for crystals smaller than 3 mm across 
is point-like, the problems associated with demagnetiza- 
tion factors is avoided. Furthermore, the magnetization 
was found to be independent of orientation. 

In order to determine accurately the critical tempera- 
ture of the ferromagnetic transition, and the critical 
exponent 7 which defines the temperature dependence 
of the zero-field magnetic susceptibility x in the critical 
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region just above the transition x(T) ^ {T/Tc — 1)"''', 
we have measured the static magnetic susceptibihty in a 
temperature interval 16 - 17.6 K. x was determined as 
a slope of the magnetization versus field curve M{H), 
through ten equidistant points between -5 and 5 Oe ev- 
ery 50 mK. To determine Tc and 7 from the experi- 
mental results we plot the inverse logarithmic derivative 
(dlnx/dT)"^ ~ —{T/Tc — versus reduced temper- 
ature e = (T — Tc)/Tc. By varying the data are 
made to appear on a straight line pointing to the ori- 
gin. The best fit of the data is shown in Fig. la for 
Tc = 16.05. From the slope we determine the exponent 
as 7 = 1.22 ± 0.02, the error reflecting variations in Tc 
within the range 16.05-16.15 where the accuracy of fit 
can not be further improved. 

0.004 



0.003 



1 0.002 



0.001 



0.000 





(b) 






P=0.75 ^ 




y=0.0026(16.1-T)°" 









14.5 14.7 14.9 15.1 



15.5 15.7 15.9 16.1 



0.00 



-0.02 



-0.04 



-0.06 



-0.08 








0.02 



0.06 



0.08 



0.04 

FIG. 1. (a)Inverse logarithmic derivative of the static sus- 
ceptibility of TDAE-Ceo, {d{\nx)/dT)~^ ~ ~{T/Tc - l)/7 
plotted against reduced temperature (T — Tc)/Tc in the criti- 
cal region above Tc. (b) Temperature dependence of the mag- 
netization of TDAE-Ceo measured in a static magnetic field of 
20 Oe (open symbols) for temperatures below Tc. Also shown 
are values of the magnetization measured at a field of 1 Oe 
(crosses) scaled by a factor of 4.75 and extrapolated values 
for zero field (bullets) scaled by factor of 6 along vertical axis. 
Solid line: the curve fit to 1/ = B{Tc - tY, as indicated. 



The critical exponent /3 describes the temperature 
dependence of the spontaneous magnetization, Ms, in 
the critical region T ^ Tc through the relation Mg ~ 
B{1 - TjTcf. To determine /3 we have measured the 
magnetization at temperatures between 14.4 K and 16 K 
in a low magnetic field. The data are shown in Fig. lb 
for 20 Oe. In fact, the functional form of the tempera- 



ture dependence of the magnetization at low fields and 
the temperature range 15.6-16.1 K does not depend on 
the value of the field. In this temperature range and for 
small H, different isotherms appear to be straight lines 
on a logif-logM plot (see Fig. 2). This implies that dif- 
ferent M{T) curves can be scaled to a unique functional 
form as Fig. lb shows. The scaled values of magnetiza- 
tion vs. temperature for lowest measured field H=l Oe 
are seen to coincide with the ones measured at 20 Oe 
near Tc- In addition, the extrapolations of the power-law 
isotherms in Fig. 2 to the H = axis leads to the same 
T-dependence of the spontaneous magnetization near the 
critical temperature, as is suggested by the susceptibil- 
ity measurements. (The corresponding curve scaled by a 
factor of 6 is also shown in Fig. 2b.) Fitting the magneti- 
zation curve to the expected power-law behavior, we find 
the best fit with Tc = 16.1 of /? = 0.75 within statistical 
error bars of 0.03 (cf. Fig. lb). 

At T = Tc the field dependence of magnetization fol- 
lows the critical isotherm H ~ M*, where S is the ex- 
ponent for the critical isotherm. To determine the criti- 
cal exponent S we have measured the magnetization over 
four decades of magnetic fields between 1 Oe and 10 kOe 
at several temperatures above and below Tc, as shown 
in Fig. 2. By fitting several isotherms near the Tc val- 
ues suggested above, we find the exponent S in the range 
S = 2.14 — 2.41, the uncertainty reflecting the chosen 
value of Tc, giving 6 = 2.28 ± 0.14. 

An important feature of the measurements is that the 
critical exponents 7, and S do not obey the scaling 
relation 7 = P(6 — 1), which is expected to apply at a 
second-order phase transition for a non-disordered sys- 
tem in equilibrium. Before we will discuss the possible 
origin of the scaling violation, we show - as an indepen- 
dent test of the consistency of the values for the critical 
exponents - that the measured data obey a general seal- 
ing form 



M{T,H) - H^^Hl {e/H^/'^\ 1 



(1) 



in the critical region at low fields H ^ Q and small rel- 
ative temperatures e ee (T - Tc)/Tc 0. Eq. @) fol- 
lows directly from the statement that the singular part of 
a thermodynamic potential (or, equivalently, its deriva- 
tives) is a generalized homogeneous function of its argu- 
ments. That is, M{b^^e,b^"H) = b^M{e,H), by taking 
b^H f{ ^ I and the standard identification of the scaling 
exponents in terms of /? and 6 (see 13|). In Fig. 3 
we plot M/H^/^ vs. X = e/H^^^^ using the values of 
the critical exponents determined above and Tc = 16.1. 
The consistency in the exponents is demonstrated by a 
"parallel" collapse of the curves for different H values. 
As discussed in Ref. the characteristic scaling func- 
tion M / H^/^ = m{x) with respect to x alone can not be 
determined directly by this fit, since neither argument of 
M(x,l) on the right-hand side of Eq. (|l|) is small in the 
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critical region. Following the procedure discussed in Ref. 

we determine the characteristic scaling function by 
plotting the reduced data m{x)/m{0) vs. x/xq, where 
xq are the values of the argument where deviations from 
power-law behavior start to occur (related to the ampli- 
tude B in the magnetization vs. temperature curve (see 
Fig. lb). The resulting plot of the data is shown in Fig. 
3b. The scaling plot confirms both the consistency of 
the measured critical exponents within the quoted error 
bars, and determines the scaling function of the phase 
transition. 
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FIG. 2. Field dependence of magnetization for different 
temperatures in the range from 14.4 K to 17.6 K (plotted 
right to left) every 0.05 K. Emphasized are several isotherms 
near Tc, in particular for T =16.1 K, 16.05 K, 16 K, 15.95 
K and 15.9 K with the corresponding fit lines for low field 
values also shown, their slopes giving the critical exponent S. 
The axes are cut off to reflect the uncertainty in the measure- 
ments. Inset: Magnetization vs. temperature data for field 
values H — 2, 3, 5, 7, 10, and 13 Oe. 



Now let us comment on the obtained values for the crit- 
ical exponents at the ferromagnetic transition in TDAE- 
Ceo (a complete list of exponents is given in Table 1). 
From previous measurements of single crystals it ap- 
pears that TDAE-Ceo for 10K< T <16K is a ferromag- 
net with localized magnetic moments with a very low 
anisotropy. Therefore, we first discuss our results in com- 
parison with the isotropic three-dimensional Heisenberg 
model. The theoretical values calculated with the Renor- 
malization Group (RG) techniques for Heisenberg 
model with spins S = 1/2 and nearest- neighbor inter- 
actions are also shown in Table 1. It is clear that the 
measured exponents differ significantly from the ones of 
3d Heisenberg model. In addition, violation of the scaling 
relation, i.e. 7 ^ l3{6 — 1), indicates an entirely different 



nature of the transition in TDAE-Cgo- In particular, an 
additional exponent 7 = (3(6 — 1) can be defined, which in 
turn violates the hyperscaling relation JlSf by an amount 
57. A modified hyperscaling relation then holds: 



2/3 + 7 = id + n)iy 



(2) 



where d = 3 is spatial dimension of the system and i> 
is the correlation length exponent. Physical insight of 
the relation (^ can be achieved by considering another 
exponent f3 defined by pi] 



(3) 



such that, together with 7/1/ = j/i/ + il the original HS 
relation is satisfied, i.e., 2/3 + 7 = di^. ^From the known 
7=0.96 we find Qv = -0.26 and thus (3 = 1.01. This 
immediately gives v — (2/3 -I- 7)/3 = 0.99. Therefore, 
51 = —0.26 < 0, meaning that the effective dimension 
d' = d + U in Eq. (^) is reduced by ~ 1/4, within the 
error bars of the measurements. A reduced effective di- 
mension indicates that the fluctuations at the transition 
are stronger than purely thermal fluctuations. For in- 
stance an enhancement of fluctuations due to conflgura- 
tional disorder in random fleld systems such as consid- 
ered in ref. For instance, within the random-fleld 
Ising model, the dimensional reduction for d — 3 was es- 
timated to lie between 57 = —1 and 57 = —1.5 |l^. The 
origin of disorder in TDAE-Cgo can indeed be related 
to random spatial realizations of the contact conflgura- 
tions of Ceo molecular orientations. As discussed in Ref. 

only the conflguration with alternating /-// orien- 
tations along c-axis is compatible with ferromagnetism, 
but other contact conflgurations may also occur with a 
flnite probability. Therefore, it is quite plausible that 
the long-range ferromagnetic order sets-in when a per- 
colating cluster of the "right" contact configurations is 
established. Indeed, the value of the exponent /3 « 1 is 
quite compatible with the backbone percolation jl^ of a 
random incipient cluster. 

Beside the measured exponents /3, 7 and S in Table 1 
we have computed the rest of the exponents using the 
valid scaling relations. Of course direct measurements of 
the exponents v, 7, rj, e.g., by scattering experiments, 
and measurements of the specific heat exponent a in 
TDAE-Cgo are necessary in order to independently con- 
firm the accuracy of these values. Preliminary measure- 
ments in fact suggest that no anomaly occurs in the spe- 
cific heat [|l^, in agreement with the predicted a < 0. 
Also, for comparison, in Table 1 we have quoted the ac- 
cepted values of the exponents for random-field Ising sys- 
tems as well as for the case of a random-exchange (or di- 
luted) Ising model ||2^ , in which disorder of the kind ex- 
pected here is a relevant perturbation at the phase tran- 
sition. At this point we wish to emphasize the distinc- 
tion between the microscopic origin of the FM interaction 
between Ceo molecules and the fluctuation mechanisms 
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leading to collective ferromagnetic behavior. The global 
response studied here does not allow us to make detailed 
conclusions about the microscopic picture, for which local 
microscopic probes such as ESR |^ and X-ray structure 
M are more appropriate. 
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FIG. 3. (a) Scaling collapse according to Eq. (1) of the 

magnetization vs. temperature curves for several values of 
the field H = 2, 3, 5, 7, 10, and 13 Oe, which are shown in 
the inset to Fig. 2. e = {T-Tc)/Tc with Tc = 16.1 and values 
of the exponents /? — 0.75 and 5 — 2.28 are used in the fit. 
The deviations from the master curve behavior for different H 
values indicate a crossover value a;o of a; = e/H^^^^ where the 
two variables in Eq. (1) become comparable, (b) The charac- 
teristic scaling function m(x)/m(0) vs. x/xq computed from 
data shown in Fig. 3a as a function of the variable x. 



TABLE I. ° Numerical values of the critical exponents of 
TDAE-Ceo: measured in this work 7, /3, and 5, and the re- 
maining exponents are computed via valid scaling relations, 
as explained in the text. Q. represents effective dimensional 
reduction in the hyperscaling relation. ^ Exponents of ran- 
dom percolation, Ref. [17]. Here Q. and /3 are the exponents 
related to the backbone percolation. Critical exponents for 
the random-field Ising model from first reference in [14]. 
Exponents for the random-exchange Ising model from Ref. 
[20], and for pure Heisenberg model, from Ref. [13], all for 
spatial dimension d = 3. 



Sys. 


7 


/3 


5 


V 


Of 




/3 


TDAE" 


1.22 


0.75 


2.28 


1.06 


-0.72 


-0.26 


1.01 


RP" 


1.82 


0.41 


5.43 


0.88 




-0.66 


0.99 


RFIM'= 


1.9 


0.06 


? 


1.02 


-0.02 


-1 


1.06 


REIM" 


1.34 


0.35 


4.78 


0.68 


-0.05 







PHM'= 


1.38 


0.36 


4.80 


0.70 


-0.12 








To conclude, the excellent reproducibility of the mea- 
sured exponents in different crystals and over time 
strongly suggests that the exponents are intrinsic to the 
material. Violation of the hyperscaling relation that im- 
plies an effective dimensional reduction by approximately 
a quarter {d! « d — 1/4) indicates that additional degrees 
of freedom — rotation of the Cgo molecules — significantly 
alters the nature of the ferromagnetic phase transition 
in TDAE-Cgo- The transition appears to be in a new 
universality class which shares some similarity with a 



backbone percolation and is attributed to the presence 
of disorder in the Cgo molecular orientations near the 
transition temperature. 
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